1 

Subwavelength Nanopatch Cavities for 
Semiconductor Plasmon Lasers 

Christina Manolatou, Farhan Rana 
School of Electrical and Computer Engineering, 
Cornell University, Ithaca, NY 14853 

Abstract 

We propose and analyze a family of nanoscale cavities for electrically-pumped surface-emitting semiconductor 
lasers that use surface plasmons to provide optical mode confinement in cavities which have dimensions in the 100-300 
nm range. The proposed laser cavities are in many ways nanoscale optical versions of micropatch antennas that are 
commonly used at microwave/RF frequencies. Surface plasmons are not only used for mode confinement but also for 
output beam shaping to realize single-lobe far-field radiation patterns with narrow beam waists from subwavelength 
size cavities. We identify the cavity modes with the largest quality factors and modal gain, and show that in the 
near-IR wavelength range (1.0-1.6 ^m) cavity losses (including surface plasmon losses) can be compensated by the 
strong mode confinement in the gain region provided by the surface plasmons themselves and the required material 
threshold gain values can be smaller than 700 cm~^. 

Index Terms 

Semiconductor Lasers, Integrated Optoelectronics, Nanotechnology, Plasmons 

I. Introduction 

Electrically pumped semiconductor lasers with nanometer scale optical cavities could be important for applications 
that benefit from ultrasmall coherent light sources, such as on-chip optical interconnects, dense photonic VLSI 
circuits, and biological or chemical sensors for micro- and nano-systems. Two questions that are interesting in 
this context are: (a) what are the smallest achievable dimensions of an electrically pumped semiconductor laser 
consistent with the current material and fabrication constrains, and (b) what is the quality of output beams shapes 
obtainable from subwavelength laser cavities. 

In the past few years much progress has been made in tightly confining light in high quality factor optical 
microcavities and in defects in ID and 2D photonic crystals 111, 12, lEI- Modal volumes close to the diffraction 
limit of (A/2n)'^ (where A is the mode wavelength and n is the refractive index seen by the mode) have been 
achieved in some of these structures ||T|, 0, Q. Optically and electrically pumped photonic crystal defect lasers 
with modal volumes few times the diffraction limit have also been demonstrated H, Q, ijS). Feedback structures, 
such as Bragg reflectors and photonic crystals, needed to achieve such small modal volumes made the overall 
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size of the laser structure several times larger than the wavelength since at least a few periods of the feedback 
structure were required for adequate mode confinement |[T1, ||2l, O, ||4l, ||5l, |l6l. Plasmonic structures for photonic 
applications have been extensively studied in the last few years lISl. Il9l. ifTOll. ifTTI. ifTH. lfT3l. lfT4l. ifBll. lfT6l. ifTTl. 
ifTsl . The large wavevector values of plasmon-polaritons near the surface plasmon resonance frequency have been 
used to achieve subwavelength device dimensions. Surface plasmon confined optical modes in waveguides have 
exhibited modal loss values ranging from 0.3 dB//im to 30 dB/fim in the visible-to-near-IR wavelength range lfT4l . 
IfTsl . While large wavevector values, and small device sizes (compared to the free-space wavelength), are possible 
for frequencies near the surface plasmon resonance frequency (which corresponds to wavelengths in the 0.4-0.6 
/im range for most important metals, such as Silver and Gold), the losses are also higher at frequencies close 
to the surface plasmon resonance frequency ||9], ifTOl - In the near-IR 1.0-1.6 /im wavelength range, although the 
losses are smaller, the wavevector values of plasmon-polaritons are also smaller ||9l, lITOl . Plasmon propagation in 
waveguides coupled to gain media has also been studied theoretically and values between 500 cm^^ and 5000 cm^^ 
for the material gain required for lossless propagation have been reported llT3l . |fT9l . An advantage of operating at 
frequencies much smaller than the surface plasmon resonance frequency is that the plasmon fields are not strongly 
confined near the surface of the metal and can therefore have significant overlap with an external gain medium. For 
realizing lasing in surface plasmon confined nanoscale optical cavities, the gain needs to not only compensate for 
intrinsic cavity losses (including surface plasmon losses) but also for losses due to external radiation. In addition, 
the output radiation patterns need to be well behaved for practical applications. In the mid-IR and far-IR wavelength 
range, where surface plasmon losses are considerably smaller compared to those at visible and near-IR wavelengths, 
surface plasmon mode confinement in dual-metal waveguides has been used to achieve lasing in quantum cascade 
devices ll20l . ||2TI . Il22l . However, these laser structures did not have subwavelength sizes in all three dimensions. 

In this paper, we propose and analyze a family of nanoscale optical cavities for electrically-pumped surface- 
emitting near-IR semiconductor lasers - semiconductor nanopatch lasers (or SNLs) - that have dimensions in the 
few hundred nanometer range and cavity volumes (not just modal volumes) approaching {X/2n)^ ||231 . Surface 
plasmons are used not only for mode confinement but also for output beam shaping to realize single-lobe far- 
field output radiation patterns with narrow beam waists. The lasers discussed here are in many ways nanoscale 
optical versions of micropatch antennas (or microstrip patch antennas) that are commonly used at microwave/RF 
frequencies ll24l . We show that in the near-IR wavelength region (1.0-1.6 /im) cavity losses in nanopatch lasers 
(including surface plasmon losses) can be compensated by gain from conventional III-V materials because surface 
plasmons can themselves be used to provide large overlap of the cavity mode with the gain region. The material 
threshold gain values needed to achieve lasing can be smaller than 700 cm^^. Despite surface plasmon losses, the 
external radiation efficiencies are between 10%-30%. Compared to all-dielectric microcavity semiconductor lasers 
reported in the literature E), ||5l, ||6l, 0, the proposed nanopatch lasers have much smaller cavity volumes and 
well behaved output beam shapes. 

In Section |ll] we discuss the basic structure of nanopatch laser cavities. Section |lll] describes the general char- 
acteristics of the modes most favorable for achieving lasing in subwavelength nanopatch cavities. Section |IV] is 
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Fig. 1. Examples of nanopatch laser structures: (LEFT) a circular nanopatch laser, (MIDDLE) a rectangular nanopatch laser, and (RIGHT) a 
hexagonal nanopatch laser. The structures are not drawn to scale. 

the main section of this paper and presents simulation resuhs for cavity modes and corresponding radiation fields. 
The focus is mainly on circular nanopatch cavities as they are the easiest to analyze computationally and the main 
characteristics of other nanopatch laser cavities are similar. Simulation results for the threshold material gain and 
threshold current values required for achieving lasing are presented in Section |V] Finally, the challenges related to 
the fabrication of nanopatch laser cavities are discussed in Section IVII 

II. Nanopatch Laser Structures 

Few examples of nanopatch laser structures are shown in Fig. [T] The basic nanopatch laser structure consists 
of a bulk semiconductor gain medium in the form of a p-i-n heterostructure sandwiched on both sides by metal 
layers that confine the lasing optical mode via surface plasmons and also serve as electrical contacts for current 
injection. The light is radiated out from the sides of the cavity. The bottom metal layer acts like an antenna reflector 
and directs the light radiated out from the sides of the cavity in the upward direction thereby contributing to the 
surface emission characteristics of nanopatch lasers. Ground planes in micropatch antennas used at microwave/RF 
frequencies perform similar functions II24I . The radii of the nanopatch lasers are between 100 nm and 300 nm and 
the heights of the dielectric part of the cavity are between 100 nm and 250 nm. The metal layers are assumed to 
be thick enough to prohibit light transmission. In the next Section, we discuss the optical modes most suitable for 
achieving lasing in nanopatch cavities. 

III. Optical Modes in Nanopatch Cavities: Discussion 

The focus in this paper will be only on those modes supported by nanopatch cavities that exhibit the smallest 
surface plasmon and radiation losses for subwavelength cavity dimensions, require the smallest cavity dimensions 
for a given lasing frequency, and have the largest overlaps with the gain region. Numerical simulations (to be 
discussed below) show that in all the nanopatch cavities discussed in this paper the most favorable modes with 
respect to the above figures of merit are similar in shape and profile. Below, we discuss these modes in more detail 
in the context of circular nanopatch lasers. 
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Fig. 2. Fields lines for (a) TMoii, (b) TMm, and (c) TMno modes of a circular nanopatch cavity. The figures on the top show the side 
view for the fields lines and the figures on the bottom show the top view in the z = h/2 plane. 



A. Optical Modes in Circular Nanopatch Cavities 

Circular nanopatch cavities are the easiest to analyze both computationally and analytically. The continuous 
angular symmetry implies that the angular dependence of the modal field amplitude can be written as exp(±ir7i(/)) 
where m is an integer For m 7^ there are two degenerate modes for each value of m corresponding to positive and 
negative values of the sign in the exponential. Equivalently, the angular dependence of the modes can be assumed 
to be cos(to0) or sin(r7i(/)) instead of exp(±im0). With reference to the vertical direction, the modes can be further 
classified as quasi-TM or quasi-TE 1241 . The TE modes do not couple well with the surface plasmons and require 
larger cavity dimensions for the same mode wavelength compared to the TM modes. Therefore, TE modes will not 
be considered in this paper The TM modes of circular nanopatch cavities bear a close resemblance to the TMj„„p 
modes of circular micropatch antennas which inside the cavity can be written as 1241 . 

E^{r, (j), z) = Eo JraiKr) cos y-j^zj eyip{±im(j)) (1) 

Er{r, (j), z) ^ Eo J'^ikr-r) sin (^^) exp(±im(/)) (2) 

/ J. \ ■ V^/^T^ JmiKr) . /pn \ , , . ,s 
E^{r, (p,z) = ^im — — Eo — sm I —z I exp(±iTO0) (3) 

where, h is the height of the cavity, fc,^ = (cj„i„p/c)^e — {pn/Kf', ujmnp is the cavity resonance frequency, and e is 
the dielectric constant of the medium inside the cavity. Although the above expressions are not exact for the TM 
modes of circular nanopatch cavities, they capture the approximate profiles and symmetries of the actual nanopatch 
cavity modes at near-IR wavelengths that are not close to the surface plasmon frequency. The TMno mode is the 
mode of choice in circular micropatch antennas at RE/microwave frequencies. This mode has a constant electric 
field component in the vertical direction and a very small radiation-Q. The small radiation-Q makes it suitable 
for antenna applications but not for achieving lasing. The electric field lines for the TMon and TMm modes of 
a nanopatch cavity are also depicted in Fig. |2] For these modes, the horizontal component of the electric field is 
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Fig. 3. (LEFT) The experimental data (SOLID WITH CIRCLES) of Ref. | 26| for the real and imaginary parts of the dielectric constant of 
Silver and the fit (DASHED) using the model plasma dispersion in the 1.0-1.8 fim wavelength range are shown. (RIGHT) The propagation 
loss of free plasmons on Silver, and the plasmon field penetration depth in the metal, are shown as a function of the wavelength. The metal is 
assumed to be surrounded by a material of refractive index 3.3 corresponding to the contact/cladding material (InGaAsP) in nanopatch lasers. 



symmetric with respect to a horizontal plane passing through the middle of the cavity and the vertical component 
of the electric field is anti-symmetric with respect to the same plane. These properties of the modes facilitate large 
mode overlap with the gain region and also result in reduced radiation losses. The reduced radiation losses are 
due to the partial cancellation of the outgoing radiation as a result of the up-down anti-symmetry of the vertical 
component of the electric field at the edges of the cavity. For a given frequency, the TMm mode results in smaller 
cavity dimensions and larger overlap of the mode with the gain region. Higher order TM modes, and in particular 
TM modes with m > 1 (whispering gallery modes), do exhibit reduced radiation losses but also require larger 
cavity dimensions and will therefore not be considered in this work. TMm mode is therefore the most suitable 
mode for achieving lasing in subwavelength circular nanopatch cavities. As discussed later, the modes most suitable 
for achieving lasing in other nanopatch cavities with subwavelength dimensions resemble the TMm mode of a 
circular nanopatch cavity. In the following Section, we present simulation results for the TMm-like cavity modes 
of nanopatch lasers. 

IV. Optical Modes in Nanopatch Cavities: Simulation Results 

A. Simulation Methods and Techniques 

For numerical simulations a 3D full-vector complex eigenmode solver, with reflectionless perfectly matched layer 
(PML) boundary conditions for outgoing radiation, is used. A 3D full-vector finite difference time domain (FDTD) 
simulator with PML boundary conditions is also used ||25| . A non-uniform computational mesh is employed and the 
minimum separation between adjacent mesh points in the metal is 4 nm, which is much smaller than the expected 
penetration depth of the plasmon field in the metal layers (see Fig. |3]l. The refractive index of the semiconductor 
gain region is assumed to be 3.5 and corresponds to the (approximately wavelength independent) refractive index of 
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InGaAsP gain region in the 1.1-1.6 wavelength range ||27l . The refractive index of the contact/cladding layers 
is assumed to be 3.3 and also wavelength independent (in the 1.0-1.6 /im wavelength range a refractive index value 
of 3.3 can always be obtained for an InGaAsP layer by adjusting the composition lIZTl ). The material loss due to 
free carriers is assumed to be 30 cm^^ per 10^* l/cm^ carrier density for holes and 5 cm^^ per 10^^ l/cm'^ carrier 
density for electrons ll28l . The n- and p-doped contact/cladding layers on each side of the gain region are both 
assumed to have carrier densities of 4 x 10^* 1/cm^. 

In all the simulations, the metal is assumed to be Silver and the experimentally measured wavelength dependent 
plasma dispersion of Silver reported by Johnson and Christy ||26l is incorporated in the simulations. In the FDTD 
method, the plasma dispersion is included by fitting the experimentally measured dispersion with the model plasma 
dispersion relation, 

, ^ 

e{uo) = - , l-, , (4) 
a; (w + z/t) 

and then adding to the FDTD equations an additional equation for the material polarization current density according 
to the model plasma dispersion, as described in Ref. ||25l . The values of Coo, ^p, and r that best fit the experimental 
data in the 1.0-1.8 /im wavelength range are 1.0, 1.38 x 10^^ rad/s, and 33 fs, respectively. The quality of the fit is 
shown in Fig. [3] Fig. [3] also shows the propagation loss of free plasmons on Silver, and the plasmon field penetration 
depth in the metal, as a function of wavelength and obtained using the model plasmon dispersion relation. These 
results will be used in the discussion that follows. In the case of the eigenmode solver, the eigenequation that needs 
to be solved numerically is. 



V X ,^ , V X H{r) = — H(r) (5) 



e(r,uj) c 

Since the dielectric constant also depends on the mode frequency, an iterative scheme is implemented in which the 
value of the dielectric constant is updated in every iteration until convergence is obtained. With PML boundary 
conditions, the eigenmode solver gives a complex mode frequency where the imaginary part is related to the photon 
lifetime in the cavity. In the FDTD method, the photon lifetime is obtained directly by observing the decay of the 
field inside the cavity with time. The photon lifetime Tp is related to the cavity quality factor Q by, Q ~ ujTp, 
where uj is the mode frequency. The photon lifetime Tp can be written as, 

111 

Tp 

where, Ti and Tr represent the intrinsic cavity loss and the radiation loss, respectively. The value of Tr is found by 
calculating the total power Pr radiated by the cavity (in the far-field) and dividing it by the total energy We stored 
in the cavity, 

1 _ Pr 

Vr ~ Wl 



(7) 



(6) 

Tr) Ti Tr 
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Fig. 4. The computed 2:-component(top view) (LEFT), z-component(side view) (MIDDLE), and the a;-component(side view) of the electric 
field for the TMm mode of a circular nanopatch cavity of height 220 nm and radius 160 nm are shown. The z-component (top view) is shown 
in a horizontal plane just below the top metal. The z- and the x-components (side view) are shown in a vertical plane passing through the 
center of the cavity. The contrast has been reduced to show the weak fringing fields in the air region. 



It needs to be pointed out here that a sHghtly different expression for the electromagnetic energy density in a 
dispersive and lossy material is given in Ref. IfTTI . However, for the parameter values chosen for the dielectric 
dispersion the difference is negligible. The integral for the total power is performed over a hemispherical surface 
of radius sufficiently large that the result is not affected by the near-fields of the cavity. The volume integral for the 
energy stored in the cavity is performed over a region slightly larger than just the dielectric and the metal regions 
of the cavity to include the fringing fields in the air region as well. We have found that in all the cases considered 
in this paper more than 90% of the field energy is stored in the dielectric and the metal regions of the cavity. Once 
Tr is found, Ti is obtained using Equation (|6]l. The external radiation efficiency Tye^jt of the laser is defined as. 

Radiation Loss Ti 

Vext = rr ■ 1 T ~ 

iotal Loss Ti + Tr 

Finally, the fraction Fg of the cavity energy in the gain region is calculated as. 



Real 



/ / I,a^n (e + {u/2)de/du;) E{r)-E*{r) dV 



Real 



(9) 



/ / La.uy {e + {ujl2)de/du:) E{r).E*{r) dV 
The results obtained from the FDTD and the eigensolution techniques were compared for several different cavities 
and the maximum difference in the calculated mode frequencies and photon lifetimes was found to be less than 
2%. The numerical error in the simulation results is expected to be of the same order 

B. Simulation Results for Circular Nanopatch Cavities 

In this Section we present simulation results for the TMm modes of circular nanopatch cavities. In the simulations 
the cavity height h and the radius R are varied. The thickness of the top and bottom contact/cladding layers is kept 
fixed and equal to 30 nm each. The contact/cladding layers need to be thick enough to incorporate the depletion 
regions (without getting fully depleted) at the interface with the smaller bandgap gain medium and also the dipole 
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Fig. 5. (LEFT) The computed wavelengths for the TMm mode of a circular nanopatch cavity are plotted as a function of the cavity radius 
for different cavity heights. The SOLID lines are the results from numerical simulations and the DASHED lines correspond to the analytical 
expression in Eq. UOt . (RIGHT) The fraction Te of the cavity electromagnetic energy in the gain region for the TMm mode is plotted as a 
function of the cavity radius for different cavity heights. The cavity heights are 140, 170, 200, 220, and 240 nm in each case. 
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Fig. 6. The computed cavity photon lifetimes (LEFT), and the fraction Fe of the cavity electromagnetic energy in the gain region (RIGHT), 
for TMm modes of circular nanopatch cavities are plotted as a function of the cavity radius for different cavity heights. The cavity heights 
are 140, 170, 200, 220, and 240 nm. The corresponding mode wavelengths are given in Fig. |5] 



layers at the interface with the metal layers. The assumed large doping level of 4 x 10^* l/cm^ for both the n- and 
p-type contact/cladding layers is therefore necessary. The thickness of the gain region is the cavity height minus 
the thickness of the two contact/cladding layers. The thickness of the top metal layer is assumed to be 100 nm. 
The computed fields for the TMm mode are shown in Fig. H] and show the large confinement of the horizontal 
component of the field in the gain region. Fig. |5] shows the wavelengths of the TMm modes as a function of the 
cavity radius for different cavity heights. Modes with wavelengths in the entke 1.0 /im to 1.6 iim range are possible 



9 




Fig. 7. The computed external radiation elSciency rjext (LEFT), and tlie cavity volume (in units of A^) (RIGHT), for the TMm modes of 
circular nanopatch cavities are plotted as a function of the cavity radius for different cavity heights. The cavity heights are 140, 170, 200, 220, 
and 240 nm. The corresponding mode wavelengths are given in Fig. \5\ The cavity volume approaches close to the diffraction limit for small 
cavity radii. 



with the cavity dimensions presented in Fig. |5] In the 1.0-1.8 fim wavelength range the interaction between the 
surface plasmons and the confined optical mode is weak and the mode shape does not deviate significantly from the 
one given by Eqs.([T]l-(|3]l- Consequently, a semi-analytical expression for the wavelength of the TM,„„p modes can 
be obtained using the approximate boundary condition that the radial component of the electric field must vanish 
at the periphery of the cavity (i.e. at r = R). If the n-th root of the derivative of the m-th Bessel function is £,mn, 
then the wavelength Xmnp of the TM„i„p mode can be written approximately as, 

_ 27rye 

Here, e is the average dielectric constant of the cavity and A is the penetration depth of the field in the top and 
bottom metal layers. Fig. |3] indicates that the penetration depth is in the 19-21 nm range for free plasmons in 
the 1.0-1.8 /im wavelength range. Simulations confirm these numbers in the case of nanopatch cavities as well. 
Assuming a constant wavelength independent penetration depth of 20 nm, an average cavity index of 3.45, the 
analytically calculated mode wavelengths (DASHED lines) for the TMm modes are compared with the numerical 
simulation results (SOLID lines) in Fig. |5] The agreement is good except for small radii when the energy in the 
fringing fields is not small. Fig. |5] also displays the fraction Fg of the cavity electromagnetic energy confined in 
the gain region as a function of the cavity radius for different cavity heights and shows that between 55%-85% of 
the mode energy is confined in the gain region. 

In plotting the simulated characteristics of the cavity modes one can choose to plot them as a function of the 
cavity radius or the mode wavelength since there is one-to-one correspondence between the two, as shown in Fig.|5] 
In the plots discussed below, we have chosen the former scheme. The corresponding mode wavelengths can be 
referenced from Fig. |5] Fig. |6] shows the cavity photon lifetimes Tp and the cavity quality factors Q for the TMm 
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mode. The resonances visible in Fig|6] appear when the diameter of the cavity is sHghtly smaller than A/2 (A is 
the mode wavelength). Fig. |7] shows the external radiation efficiency r/ext- For small cavity radii and larger cavity 
heights the radiative losses dominate, and for large cavity radii the intrinsic losses dominate. In order to obtain 
decent external efficiencies the radiative losses cannot be much smaller than the intrinsic losses. Therefore, there is 
little to gained in terms of laser performance by reducing the radiative losses much below the intrinsic losses. An 
interesting aspect of the nanopatch cavities is their small cavity size. Fig. |7] shows the cavity volume (not the mode 
volume) in units of A'^ as a function of the cavity radius. For comparison, the diffraction limit for the mode volume 
{X/2n)^ is also indicated in the figure (using an average index value of 3.45 for the dielectric part of the cavity). 
Although the smallest achievable cavity volume using plasmon confinement is neither related to, nor constrained 
by, the diffraction limit {X/2n)^, the diffraction limit has become a useful standard for measuring Ught confinement 
in small cavities ||291 . Fig. |7] shows that the cavity volume approaches values close to the diffraction limit for small 
radii. 



C. Far-Field Radiation Patterns of Circular Nanopatch Cavities 

A useful property of the proposed nanopatch lasers is their surface-normal emission characteristics. Nanopatch 
lasers can emit in the surface normal direction with a single-lobe far-field radiation pattern and a relatively narrow 
beam waist. The radiation pattern p{9, (p) is defined as, 

Pie,ci>)^^^^^iMl (11) 

S{r, 6, 0).r|max over 6,4, 

where S{r, 9, cj)) is the Poynting vector. In simulations, the radiation pattern is found by finding the Poynting vector 
from the calculated fields over a hemispherical surface of radius large enough such that cavity near-fields do not 
affect the results. Typically, a value of radius between 4-6 fim is an adequate compromise between accuracy and 
Umited computational resources. For a circular nanopatch cavity, the radiation pattern will not have a 0-dependence 
if the modes are assumed to vary with (p as exp(±im0). On the other hand, if the 0-dependence of the modes 
is assumed to be cos(m ip) or sin(m (p) then the radiation pattern will also exhibit a 0-dependence. The radiation 
pattern p{9) for exp(±i m cp) modes can be obtained by averaging the radiation pattern p{9, (p) for cos(to (p) or 
sm{m(p) modes over (p. Fig.[8]shows the far-field radiation pattern for the TMm mode of a 160 nm radius and 220 
nm cavity height circular nanopatch cavity. The angular dependence of the z-component of the field is assumed to 
be proportional to cos{(p) rather than exp(±i (p) to show the 0-dependence of the pattern. The radiation is emitted 
out from the sides of the cavity but interferes constructively in the upward direction. This is also shown explicitly 
in Fig. |9] A weak surface plasmon wave is also emitted along the bottom metal layer The surface wave decays as 
it propagates and its power would have made no contribution to the radiation pattern if the radiation pattern were 
computed over a hemispherical surface of much larger radius. 
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Fig. 8. A circular naiiopatch laser (LEFT) and the computed far-field radiation pattern (RIGHT) for the TMm mode of a 160 nm radius and 
220 nm cavity height laser are shown. The angular dependence of the z-component of the field is assumed to be proportional to cos((/)) rather 
than exp(ibi (p). A weak surface plasmon wave emitted along the bottom metal layer is also visible in the radiation pattern. 
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Fig. 9. The computed z-component (LEFT) and the 2-component (RIGHT) of the radiated field for the TMm mode of a circular nanopatch 
cavity of radius 160 irai and height 220 nm are shown. The angular dependence of the 2-component of the field is assumed to be proportional 
to cos(</)) rather than exp(±i 0). The radiation fields are shown in the E-plane (ip = 0)- The contrast in the figures has been greatly reduced 
to show the radiated fields compared to the much stronger fields inside the cavity. 

D. Analytical Model for the Far-Field Radiation Patterns of Nanopatch Cavities 

The radiation patterns, and the surface-normal emission of radiation in particular, can be understood and estimated 
from considerations similar to those used in the analysis of aperture antennas at RF/microwave frequencies |[24l . 
The sides of the nanopatch cavities can be considered as apertures from which radiation takes place. In order to 
calculate the radiated fields, equivalent electric and magnetic surface current densities, K{r) and M{r), respectively, 
given by, 

K{r) =nX H(r) I boundary (12) 
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M(r) = -n X £;(r)|boundary (13) 

can be assumed to exist at the surface of the cavity and the radiated fields can be assumed to be generated by these 
current densities ||24l . Here, n is a unit vector perpendicular to the surface of the cavity. Since the tangential magnetic 
fields at the sides of the nanopatch cavities are very small (and equal to zero from Eqs.dT]!-©), the electric current 
density can be ignored for simplicity. Unlike at RF/microwave frequencies, the metal layers at near-IR frequencies 
cannot be considered as perfectly conducting since the real part of the dielectric constant e{Lo) dominates at near-IR 
frequencies, as shown in Fig. |3] However, the large negative values of e(a;) at near-IR frequencies imply that the 
boundary condition of zero tangential electric field at the surface of the metal holds approximately. Therefore, the 
magnetic surface current density M{r) will have its image in the bottom metal layer. For the TMm mode of a 
circular nanopatch cavity, the components of the magnetic surface current density and their images are depicted in 
Fig. [To] The surface normal emission of radiation by nanopatch cavities can be understood as the radiation emitted 
by the equivalent magnetic surface current density and its image. In the far-field, if the electric field is represented 
by a vector potential F{f) using the expression, 

then the vector potential corresponding to the magnetic surface current density M{r) can be written as, 

F(r)= ^oeM^kr) f f da' M{r')eM-ik.?) (15) 

4 7r r J J boundary 

Using the expressions for the field components of the TM modes of a circular nanopatch cavity given in Eqs.([T]l-(l3]l, 
the far-field radiation pattern p{9, (f)) for the TMm mode can be calculated analytically and the approximate result 
is, 

pi0,(l)) oc [sin(fcft. cos(6'))2fc/i cos(6') + Scos(fc/i cos(6')/2)]^cos2((/)) 
+ [{sm{kh cos{0))2kh cos{9) + Bcos{kh cos{e)/2)) cos{9) 

-Asm{kh coa{e)) sin^ (9)]^ sin^ (cj)) (16) 

where, A — 2Tr^k/{k'^h), k — lo/c, h is the cavity height, and kr = ^u/R. In deriving the above expression, 
the (/(-dependence of the z-component of the field is assumed to be cos((/)). The actual TMm mode is not exactly 
symmetric with respect to a horizontal plane passing through the center of the cavity, as assumed in Eqs.([TJ-([3j. The 
value of the dimensionless constant B is chosen to model this asymmetry and numerical simulations show that the 
value of B is in the 0.05 to 0.2 range. The simulated and the analytically calculated (using the expression above) 
radiation patterns for the TMm mode of a circular nanopatch cavity are shown in Fig.[TOl The cavity height is 220 
nm, the cavity radius is 160 nm, and the value of B is 0.15. The agreement in both the H-plane (</> — 90-degrees) 
and E-plane (cj) = 0) is good except for values of 9 close to 90-degrees in the E-plane. The error for values of 9 
close to 90-degrees in the E-plane is due to the fact that surface plasmon waves on the bottom metal layer that 
are emitted in the E-plane are not adequately described by the analytical model based on the equivalent magnetic 
current densities. 
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Fig. 10. (LEFT) The components of the equivalent magnetic surface current density and their images are depicted for the TMm mode of a 
cii'cular nanopatch cavity. (MIDDLE and RIGHT) The simulated (SOLID) and the analytically calculated (DASHED) radiation patterns p{8, 4>) 
for the TMiii mode of a circular nanopatch cavity are shown in the H-plane (<j> = 90-degrees) and in the E-plane (0 = 0). The cavity height 
is 220 nm and the cavity radius is 160 nm. The ^-dependence of the z-component of the field is assumed to be cos{<j>). The en'or for 8 close to 
90-degrees in the E-plane is due to the fact that surface plasmon waves are not adequately described by the analytical model based on equivalent 
magnetic current densities. 




Fig. 11. The simulated far-field radiation patterns p{6) for the TMm modes of circular nanopatch cavities of height 220 nm and different 
radii are shown as a function of the angle 9. The ^-dependence of the modes is assumed to be exp(±i(/>) and therefore the radiation patterns 
depend only on the angle 9. 
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Fig. 12. The computed z-component(top view) (LEFT), 2;-component(side view) (MIDDLE), and tlie a;-component(.side view) of the electric 
field for the TMm-like mode of a square nanopatch cavity of side 280 nm and height 220 nm are shown. The 2-component (top view) is 
shown in a horizontal plane just below the top metal. The contrast has been adjusted to show the weak fringing fields. 




Fig. 13. The computed z-component(top view) (LEFT), 2:-component(side view) (MIDDLE), and the a;-component(side view) of the electric 
field for the TMm-like mode of a hexagonal nanopatch cavity of height 220 nm and side 200 nm are shown. The z-component (top view) is 
shown in a horizontal plane just below the top metal. The contrast has been adjusted to show the weak fringing fields. 

Fig, fm shows the computed far-field radiation patterns p{9) for the TMm modes of circular nanopatch cavities of 
height 220 nm for different cavity radii. The ^-dependence of the modes is assumed to be exp(±i (p) for simplicity. 
It can be seen that the radiation patterns depend sensitively on the cavity dimensions. In the surface normal direction 
the radiation is entirely due to the 0-component of the magnetic current density. For small cavity radii (R < h 
nm), the radiation is dominated by the (/i-component of the magnetic current density. For larger cavity radii {R > h 
nm), the z-component of the magnetic current density also contributes strongly to the radiation. 

E. Modes, Degeneracies, and Far-Field Radiation Patterns in Square and Hexagonal Nanopatch Cavities 

In actual fabricated circular nanopatch cavities the angular symmetry is expected to be broken unintentionally and 
the modes are expected to have cos(m (/>) or sin(m ((>) dependence rather than exp(±i m(j)) dependence. The angular 
symmetry can be intentionally broken using the square and hexagonal nanopatch cavities that also support modes 
that resemble the TMm modes of circular nanopatch cavities. These modes have a horizontal field component 
that is symmetric with respect to a horizontal plane passing through the center of the cavity and a vertical field 
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Fig. 14. The computed far-field radiation patterns for the TMm-like modes of a square (LEFT) and hexagonal (RIGHT) nanopatch cavities 
are shown. The cavity dimensions are the same as in Figs. nH and [T3I 



component that is antisymmetric with respect to the same plane. The antisymmetry of the vertical component of 
the field at the cavity edges helps in reducing the radiation losses and the horizontal field component enjoys a large 
overlap with the gain region. 

The calculated fields for a square nanopatch cavity of height 220 nm and side 280 nm and a hexagonal nanopatch 
cavity of height 220 nm and side 200 nm are shown in Figs. [12] and [13] respectively. The mode wavelength and 
the cavity-(3 for the square cavity are 1.30 /im and 196, respectively, and 1.41 /im and 242, respectively, for the 
hexagonal cavity. These cavity-(5 values are comparable to those of circular nanopatch cavities. The continuous 
angular symmetry in circular nanopatch cavities implies that the two modes with the same angular index m (for 
m ^ 0) are degenerate. The symmetries in case of square and hexagonal nanopatch cavities are discrete, and 
correspond to the and C^v point groups, respectively |[30l . The character tables of the and C^v point 
groups are given in Appendix [A] The modes for the square and hexagonal nanopatch cavities shown in Figs. [T2l 
and [T3] correspond to the two-dimensional E and Ei representations of the €4^ and Cqv point groups, respectively, 
and are therefore each doubly degenerate. Most of the properties of square and hexagonal nanopatch cavities are 
not too different from those of circular nanopatch cavities. The far-field radiation patterns for a square nanopatch 
cavity of height 220 nm and side 280 nm and a hexagonal nanopatch cavity of height 220 nm and side 200 nm 
are shown in Fig. [14] and resemble the radiation patterns of circular nanopatch cavities (see Fig. [8]). 
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TABLE I 

Laser parameter values used in calculation of threshold currents ||28l 



Parameter 


Value 


Gain coefficient go 

Transparency carrier density Ntr 

Auger Coefficient C 

Spontaneous emission factor n^p 

Surface recombination velocity fsur/ 1311 

Index of active region n 


2000 cm-^ 
1.8 X 10^** cm-3 
3 X 10^29 cm'S/s 
1.5 

10^ cm/s 
3.5 



V. Threshold Gain and Threshold Current of Nanopatch Lasers 
A. Threshold Gain 

In a semiconductor laser the material threshold gain gth, the mode energy gain confinement factor Fg, the gain 
medium refractive index n, and the photon hfetime Tp are related as 1281 , 

re-gth = - (17) 

n Tp 

where, c is the speed of light in free space. The strong mode confinement in the gain region provided by the 
surface plasmons in nanopatch cavities results in large modal gain values that are needed to compensate for the 
large cavity losses. The photon lifetimes and the mode energy gain confinement factors in conventional edge- 
emitting semiconductor lasers typically range from 1-4 ps and 0.04-0.08%, respectively f28|. The results presented 
in Section ITV-B I show that the photon lifetimes in nanopatch cavities are almost 10 times smaller than in conventional 
semiconductor lasers but the mode energy gain confinement factors are almost 10 times larger compared to their 
values in conventional semiconductor lasers. Therefore, the material threshold gain values in nanopatch lasers are 
not expected to be much different than in conventional semiconductor lasers. Although Eq. [17] can be used to find 
the material threshold gain values from the calculated values of Fg and Tp, this procedure does not take into account 
modification of the mode profile in the presence of material gain. A more accurate method to obtain the threshold 
gain is to calculate the complex mode frequency in the presence of the material gain using the techniques discussed 
in Section IIV-AI and then iteratively adjust the value of the material gain until the imaginary part of the complex 
mode frequency becomes zero (or sufficiently small). This procedure has been used this paper and the results for 
the threshold gain of circular nanopatch lasers are shown in Fig. [15] The lasing mode is assumed to be TMm. The 
threshold gain values are less than 1000 cm^^ for cavity radii larger than 175 nm and cavity heights larger than 170 
nm. In the 1.5-1.6 iim wavelength range, gain values smaller than 700 cm^^ are achievable. III-V semiconductor 
materials, such as InGaAsP and InGaAlAs, can easily provide material gain values in the 1200-1600 cm~^ range 
in the 1.2-1.6 /im wavelength range and are therefore suitable for nanopatch lasers li28l . 
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Fig. 15. (LEFT) The computed material threshold gain g^^ for the TMm mode of a circular nanopatch laser for different cavity heights 
are shown. (RIGHT) The thi'eshold cun'ent 1^^ for a circular nanopatch laser are shown as a function of the cavity radius for different cavity 
heights. The cavity heights are 140, 170, 200, and 240 nm. The top and bottom contact/cladding layers are 30 nm thick each. The values of 
various parameters used in the calculations are listed in Table. 1. 



B. Threshold Current 

In this Section, estimates of the threshold currents in semiconductor nanopatch lasers are presented. The threshold 
carrier density Nth can be estimated from the expression for the material gain as a function of the carrier density 
N I28i, 

.(iV)=.olog(|-) (18) 

The threshold current Ith can be related to the electron-hole recombination rate R{N) in the active region at 
threshold, 

% = R{Nth) (19) 
qV 

Here, q is the electron charge and V is the volume of the active region. The carrier density dependent recombination 
rate R{N) (units: cm^'^-s^^) is assumed to have dominant contributions from surface recombination, spontaneous 
emission, and Auger scattering. Surface recombination is expected to become more important as the cavity radius 
decreases and the surface area to volume ratio of the active region increases. The net recombination rate R{N) can 
be written as, 

R{N) = Rsp{N) + R 

auger 

= re-9{N)'^ + CN'' + v,urfN^ (20) 
n Vp V 

where, c is the speed of light, n is the index of the active region, Usp is the spontaneous emission factor, Vp is 

the mode volume and is approximately equal to the cavity volume for nanopatch lasers, C is the Auger coefficient, 

and Vsurf is the surface recombination velocity and equals ^ 10'^ cm/s for passivated InGaAsP active layers |3l1 . 



18 



(A/V) is the surface area to volume ratio for the active region and equals 2/R for a circular nanopatch laser of 
radius R. Given the small cavity size, the spontaneous emission is assumed to occur only in the lasing mode. Using 
the parameter values given in Table. 1, the calculated threshold currents of circular nanopatch lasers are shown 
in Fig. \T5\ as a function of the cavity radius for different cavity heights. The parameter values are assumed to 
be constant whereas in reality they will change shghtly as the composition of the gain medium is changed. The 
parameter values shown in Table. 1 correspond to bulk InGaAsP active medium lattice matched to InP and with a 
bandgap corresponding to 1.3 /im wavelength ||28l . The threshold currents are large for very small radii because 
the photon lifetimes are small. For large cavity radii the threshold currents are again large because the active region 
volumes are large. Therefore, there is an optimum value of the cavity radius that results in the minimum value of 
the threshold current. The minimum values of the threshold current are in the 3 — 5 range and correspond to 
threshold current densities in the 2-3 kA/cm^ range. 

VI. Conclusion 

Challenges associated with the fabrication of nanopatch lasers are not expected to be impossible to meet. Substrate 
removal techniques used in the fabrication of dual-metal-waveguide semiconductor far-IR lasers can be employed 
for the realization of the dual-metal nanopatch laser structures ll32l . However, in the case of nanopatch lasers the 
semiconductor films obtained after substrate removal would be much thinner and would require more care. The 
metal used for the top and bottom layers must make good ohmic contacts to the contact/cladding layers, must not 
react with or diffuse into the semiconductor, and must also have small surface plasmon losses. Silver seems to be 
the ideal choice. The work function of Silver is very close to that of Titanium and Titanium/Gold ohmic contacts 
are commonly used for both n- and p-type III-V semiconductors f33l . Good quality Silver ohmic contacts, which do 
not require post-deposition rapid thermal anneal, to moderately doped InP have been demonstrated 1341 . Electrical 
contact leads to the top metal layer would need to be microfabricated since the laser structure is too small to be 
contacted directly with electrical probes. Finally, the exposed side walls of the cavity would need to be passivated to 
reduce surface recombination. Quantum well gain media can also be used in place of bulk gain media in nanopatch 
lasers. Quantum wells can generally provide larger material gain compared to bulk for the same carrier density due 
to the reduced density of states 1281 . However, a quantum well gain medium would have smaller overlap with the 
cavity mode and the material gain would also be polarization dependent. The effect of interface roughness on the 
cavity quality factors has not been considered in this work. In nanoscale lasers, the surface area to volume ratio 
increases as the cavity dimensions shrink and therefore interface roughness will likely play an important role. 

The ideas presented in this paper show that in the 1.0-1.6 /im wavelength range surface plasmons can be used 
to confine light in nanoscale optical cavities. At these longer wavelengths since the surface plasmon fields are not 
strongly confined to the surface of the metal, they can have significant overlap with an external gain medium. In 
fact, surface plasmon confinement in the 1.0-1.6 /im wavelength range can be used to obtain mode gain confinement 
factors that are substantially larger than in conventional semiconductor lasers thereby providing large modal gain 
values necessary to achieve lasing. 
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The benefits of nanopatch laser structures are not expected to be limited to just reduced dimensions. The 
small photon lifetimes in nanopatch lasers compared to conventional semiconductor lasers could enable ultrawide 
bandwidths (> 100 GHz) for direct current modulation. The dual-metal layer structure of nanopatch lasers could 
enable them to be used as functional elements and placed on desired substrates, such as SiUcon microchips, to 
provide lasers on demand capabilities for multifunctional micro- and nano-systems. 

The authors would hke to thank Dr. Filbert BartoU for support through the National Science Foundation faculty 
CAREER award and Dr. Leda Lunardi for support through the National Science Foundation award #0636593, and 
also TLX Ughtwave Corp., Intel Inc., and Infotonics Inc. for support. 
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Appendix A 

Character Tables for the and Cq^ Point Groups 

The character tables for C4„ and Cg^ point groups are shown below. The notations for the group symmetry 
operations are according to Ref. ll30l . 



6*41, Point Group 


Civ 


E 


2C4 


C2 


2av 


2CTd 


Ai 


1 


1 


1 


1 


1 


A2 


1 


1 


1 


-1 


-1 


Bi 


1 


-1 


1 


1 


-1 


B2 


1 


-1 


1 


-1 


1 


E 


2 





-2 













Cqv 


Point Group 






Cqv 


E 


2C6 


2C3 


C2 




3(Tx 


Ai 


1 


1 


1 


1 


1 


1 


A2 


1 


1 


1 


1 


-1 


-1 


Bi 


1 


-1 


1 


-1 


1 


-1 


B2 


1 


-1 


1 


-1 


-1 


1 


El 


2 


1 


-1 


-2 








E2 


2 


-1 


-1 


2 









The modes discussed in this paper in square and hexagonal nanopatch cavities correspond to the two-dimensional 
E and Ei representations of the Cav and Cqv point groups, respectively, and are therefore each doubly degenerate. 
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List of Figures 

Examples of nanopatch laser structures: (LEFT) a circular nanopatch laser, (MIDDLE) a rectangular 
nanopatch laser, and (RIGHT) a hexagonal nanopatch laser The structures are not drawn to scale. . . [s] 
Fields lines for (a) TMon, (b) TMm, and (c) TMno modes of a circular nanopatch cavity. The figures 
on the top show the side view for the fields lines and the figures on the bottom show the top view in 

the z plane 

(LEFT) The experimental data (SOLID WITH CIRCLES) of Ref. EH for the real and imaginary parts 
of the dielectric constant of Silver and the fit (DASHED) using the model plasma dispersion in the 
L0-L8 /im wavelength range are shown. (RIGHT) The propagation loss of free plasmons on Silver, 
and the plasmon field penetration depth in the metal, are shown as a function of the wavelength. 
The metal is assumed to be surrounded by a material of refractive index 3.3 corresponding to the 

contact/cladding material (InGaAsP) in nanopatch lasers Q 

The computed z-component(top view) (LEFT), 2:-component(side view) (MIDDLE), and the a;-component(side 
view) of the electric field for the TMm mode of a circular nanopatch cavity of height 220 nm and 
radius 160 nm are shown. The z-component (top view) is shown in a horizontal plane just below the 
top metal. The z- and the x-components (side view) are shown in a vertical plane passing through the 
center of the cavity. The contrast has been reduced to show the weak fringing fields in the air region. Q 
(LEFT) The computed wavelengths for the TMm mode of a circular nanopatch cavity are plotted as a 
function of the cavity radius for different cavity heights. The SOLID lines are the results from numerical 
simulations and the DASHED lines correspond to the analytical expression in Eq.(fTO]i. (RIGHT) The 
fraction Fe of the cavity electromagnetic energy in the gain region for the TMm mode is plotted as 
a function of the cavity radius for different cavity heights. The cavity heights are 140, 170, 200, 220, 

and 240 nm in each case [s] 

The computed cavity photon lifetimes (LEFT), and the fraction of the cavity electromagnetic energy 
in the gain region (RIGHT), for TMm modes of circular nanopatch cavities are plotted as a function 
of the cavity radius for different cavity heights. The cavity heights are 140, 170, 200, 220, and 240 

nm. The corresponding mode wavelengths are given in Fig. |5] Is] 

The computed external radiation efficiency rj^^t (LEFT), and the cavity volume (in units of A^) 
(RIGHT), for the TMm modes of circular nanopatch cavities are plotted as a function of the cavity 
radius for different cavity heights. The cavity heights are 140, 170, 200, 220, and 240 nm. The 
corresponding mode wavelengths are given in Fig. |5] The cavity volume approaches close to the 

diffraction limit for small cavity radii Isl 

A circular nanopatch laser (LEFT) and the computed far-field radiation pattern (RIGHT) for the TMm 
mode of a 160 nm radius and 220 nm cavity height laser are shown. The angular dependence of the z- 
component of the field is assumed to be proportional to cos((/i) rather than exp(±i 0). A weak surface 
plasmon wave emitted along the bottom metal layer is also visible in the radiation pattern [h] 
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9 The computed x-component (LEFT) and the z-component (RIGHT) of the radiated field for the TMm 
mode of a circular nanopatch cavity of radius 160 nm and height 220 nm are shown. The angular 
dependence of the z-component of the field is assumed to be proportional to cos((/)) rather than 
exp(±? (f). The radiation fields are shown in the E-plane (0 — 0). The contrast in the figures has been 
greatly reduced to show the radiated fields compared to the much stronger fields inside the cavity. . . [h] 

10 (LEFT) The components of the equivalent magnetic surface current density and their images are 
depicted for the TMm mode of a circular nanopatch cavity. (MIDDLE and RIGHT) The simulated 
(SOLID) and the analytically calculated (DASHED) radiation patterns p{0, (f)) for the TMm mode of 
a circular nanopatch cavity are shown in the H-plane ((j> — 90-degrees) and in the E-plane ((j> — 0). 
The cavity height is 220 nm and the cavity radius is 160 nm. The ^-dependence of the z-component of 
the field is assumed to be cos{<f>). The error for 9 close to 90-degrees in the E-plane is due to the fact 
that surface plasmon waves are not adequately described by the analytical model based on equivalent 
magnetic current densities [ij] 

11 The simulated far-field radiation patterns p{9) for the TMm modes of circular nanopatch cavities of 
height 220 nm and different radii are shown as a function of the angle 9. The 0-dependence of the 
modes is assumed to be exp(±i cf)) and therefore the radiation patterns depend only on the angle 9. . [ij] 

12 The computed z-component(top view) (LEFT), 2r-component(side view) (MIDDLE), and the 2;-component(side 
view) of the electric field for the TMm-like mode of a square nanopatch cavity of side 280 nm and 
height 220 nm are shown. The z-component (top view) is shown in a horizontal plane just below the 

top metal. The contrast has been adjusted to show the weak fringing fields 1^ 

13 The computed z-component(top view) (LEFT), z-component(side view) (MIDDLE), and the a;-component(side 
view) of the electric field for the TMm-like mode of a hexagonal nanopatch cavity of height 220 nm 

and side 200 nm are shown. The z-component (top view) is shown in a horizontal plane just below 

the top metal. The contrast has been adjusted to show the weak fringing fields 1^ 

14 The computed far-field radiation patterns for the TMm-like modes of a square (LEFT) and hexagonal 
(RIGHT) nanopatch cavities are shown. The cavity dimensions are the same as in Figs. [T2] and [T3] Q 

15 (LEFT) The computed material threshold gain gth for the TMm mode of a circular nanopatch laser 
for different cavity heights are shown. (RIGHT) The threshold current Ith for a circular nanopatch 
laser are shown as a function of the cavity radius for different cavity heights. The cavity heights are 
140, 170, 200, and 240 nm. The top and bottom contact/cladding layers are 30 nm thick each. The 
values of various parameters used in the calculations are listed in Table. 1 Q 
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